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We investigate universal behavior of isolated many-body systems far from equilibrium, which 
is relevant for a wide range of applications from ultracold quantum gases to high-energy particle 
physics. The universality is based on the existence of nonthermal fixed points, which represent 
nonequilibrium attractor solutions with self-similar scaling behavior. The corresponding dynamic 
universality classes turn out to be remarkably large, encompassing both relativistic as well as non¬ 
relativistic quantum and classical systems. For the examples of nonrelativistic (Gross-Pitaevskii) 
and relativistic scalar field theory with quartic self-interactions, we demonstrate that infrared scal¬ 
ing exponents as well as scaling functions agree. We perform two independent nonperturbative 
calculations, first by using classical-statistical lattice simulation techniques and second by applying 
a vertex-resummed kinetic theory. The latter extends kinetic descriptions to the nonperturbative 
regime of overoccupied modes. Our results open new perspectives to learn from experiments with 
cold atoms aspects about the dynamics during the early stages of our universe. 


I. INTRODUCTION 
A. Universality far from equilibrium 

The notion of universality or scaling phenomena far 
from equilibrium is to a large extent unexplored, in par¬ 
ticular, in relativistic quantum held theories. Here the 
strong interest is mainly driven by advances in our un¬ 
derstanding of the early universe dynamics as well as rel¬ 
ativistic collision experiments of heavy nuclei in the lab¬ 
oratory. Experimentally, the investigation of quantum 
systems far from equilibrium is also boosted by the non¬ 
relativistic physics of ultracold quantum gases. By using 
optical or atom chip traps, they provide a flexible testbed 
with tunable interactions, symmetries and dimensional¬ 
ity, with connections to a wide variety of physical sys¬ 
tems. This may include high-energy physics, in partic¬ 
ular, if both relativistic and nonrelativistic systems be¬ 
long to the same universality class [Tj. It is the aim of 
this paper to establish such a connection for bosonic held 
theories far from equilibrium. 

In recent years, important new universality classes 
have been discovered in isolated relativistic systems far 
from equilibrium PHIZ]. The universality is based on the 
existence of nonthermal hxed points ®mm, which 
represent nonequilibrium attractor solutions with self¬ 
similar scaling behavior. The underlying physical pro¬ 
cesses are reminiscent of the stationary transport of 
conserved charges in the phenomenon of wave turbu¬ 
lence m EH- However, no external sources or sinks 
are present for isolated systems. In high-energy physics, 
collision experiments of heavy nuclei (where a ‘fireball’ 
expands against the surrounding vacuum) and the evo- 


* pineiroorioli@thphys.uni-heidelberg.de 
t K.Boguslavski@thphys.uni-heidelberg.de 
t J.Berges@thphys.uni-heidelberg.de 



FIG. 1. Schematic illustration of the occupation number dis¬ 
tribution near a nonthermal fixed point, shown as a function 
of momentum p for two subsequent times ti and £2 > £ 1 . 
The scaling exponents a and /3 characterize the self-similar 
evolution according to equation 0 - 


lution of the early universe provide important examples 
for isolated quantum systems. 

Also experimental setups employing ultracold quan¬ 
tum gases can be largely isolated, such that their dy¬ 
namics is governed by unitary time evolution. Nonther- 
mal fixed points have been investigated in this context 
for Bose gases HSII22H2SI- In these nonrelativistic sys¬ 
tems, Bose condensation occurs out-of-equilibrium as a 
consequence of an inverse particle cascade J2U G23 EOHSl] • 
Similar cascades with nonthermal Bose condensation are 
also known for relativistic scalar field theories in vari¬ 
ous spatial dimensions n n E HU M- The question 
arises, whether these similarities between the relativistic 
and nonrelativistic systems can be made quantitative. 
Clearly, there are important differences and one has to 
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specify which properties can be universal. 

To illustrate this, FIG. |T] shows a typical distribution 
function f(t, p) near a nonthermal fixed point as a func¬ 
tion of momentum p for two subsequent times t = t\ and 
f 2 > t\. We distinguish different momentum ranges by 
their occupancies /(t, p) in terms of a small (interaction 
or ‘diluteness’) parameter A < 1 as will be explained 
below. The inverse particle cascade leading to Bose con¬ 
densation occurs in the highly occupied low-momentum 
regime, where /(f,p) > 1/A (3]. This particle transport 
towards low momenta is part of a dual cascade, in which 
energy is also transferred by weak wave-turbulence to¬ 
wards higher momenta (1/A /(t,p) 1). The latter 

evolves until a high-momentum scale is reached, where 
the characteristic mode occupancy becomes comparable 
to the ‘quantum-half’ [2j. 

The different cascades exhibit approximate power-law 
behavior for mode occupancies f(t, p) in characteristic 
inertial ranges of momenta. However, in general iso¬ 
lated systems out of equilibrium cannot realize station¬ 
ary transport solutions in the absence of external driv¬ 
ing forces. Instead, the transport in isolated systems is 
described in terms of the more general notion of a self¬ 
similar evolution, where the distribution function obeys 
for isotropic systems 


/(i,p) = t a fs(£ = tP\p\) (1) 


in a given scaling regime. Here, all quantities are consid¬ 
ered to be dimensionless by use of some suitable momen¬ 
tum scale, which is specified below. 

The values of the scaling exponents a and /3, as well 
as the form of the nonthermal fixed point distribution 
fs(f) are universal. More precisely, all models in the 
same universality class can be related by a multiplica¬ 
tive rescaling of t and p. Quantities which are invari¬ 
ant under this rescaling are universal. Accordingly, all 
system-dependent aspects of the distribution are con¬ 
tained in two nonuniversal amplitudes, which we denote 
as A and B. It is convenient to define them according 
to /s(£ = B) = A with dfs(£ = B)/df = —2 A/B such 
that A characterizes the amplitude of the scaling func¬ 
tion at £ = B 1 where the occupation number receives its 
dominant contribution. 

We emphasize that the universal properties can be dif¬ 
ferent for different inertial ranges. This is indicated in 
FIG. [l] where in the direct cascade regime other scaling 
exponents a' , (3' and a different scaling function f' s than 
in the infrared are found. For instance, two theories can 
have the same universal low-momentum properties while 
they may differ significantly in another inertial range at 
higher scales. This is very similar to the classification 
of universal properties in thermal equilibrium, where one 
distinguishes for a given theory between infrared and ul¬ 
traviolet fixed points and associated scaling properties 
depending on the momentum regime. 


B. Outline of results 


In this work we compute the exponents a, (3 and the 
scaling function fg of the self-similar distribution |l]) in 
the infrared regime. We present results for nonrelativistic 
(Gross-Pitaevskii) as well as (massless) relativistic scalar 
field theory with quartic self-interactions, respectively. 
While the relativistic theory captures important aspects 
of the Higgs sector of particle physics or of inflationary 
models for early universe dynamics, the Gross-Pitaevskii 
field theory can describe a dilute Bose gas. 

A central conclusion of this paper is that the infrared 
scaling exponents and scaling functions of these theories 
belong to the same universality class, i.e. the results for a, 
f3 and the universal form of the scaling function fg in (|T|) 
agree. This is nontrivial, in particular, since the nonrel¬ 
ativistic system conserves total particle number whereas 
in the relativistic theory number-changing processes are 
possible. The agreement found in the infrared is also re¬ 
markable in view of the fact that in the inertial range of 
the direct energy cascade towards higher momenta the 
exponents a ', /3' and the scaling function f' s from the 
different theories do not agree mm- 

Since the large occupancies at low momenta lead to 
strongly nonlinear dynamics, one cannot apply standard 
perturbative kinetic theory in the infrared. (Perturba¬ 
tive approaches mmm are often used to describe the 
direct energy cascade at higher momenta]^]) We perform 
two independent nonperturbative calculations. The first 
approach employs classical-statistical lattice simulation 
techniques in sections |H| and III [6}|35]. The second, ana¬ 
lytical method applies a vertex-resunnned kinetic theory 
in sections IV and [Vj which is based on an expansion 
in the number of held components N to next-to-leading 
order m i33- The approach extends well-established 
kinetic descriptions E Hi mi to the nonperturbative 
regime of overoccupied modes. 

This is used to obtain the analytic estimate for the 
scaling exponents of the 


particle transport: a = (3 d, (3 = - , (2) 

towards low momenta, with spatial dimension d. This 
is a central analytic result of this work. In contrast to 
the previously known negative values for a and j3 from 
perturbative estimates mm , the positive values obtained 
from the vertex-resummed kinetic theory describe an in¬ 
verse particle transport with growing occupation num¬ 
ber in the infrared. The fixed relation between a and 
/3 reflects the conservation of particle number density 
n = f d d p/(2n) d f(t, p) ~ in this (nonrelativis¬ 

tic) inertial range by using the self-similarity ([l]). Most 


1 Perturbative kinetic theory has also been employed in the in¬ 
frared to describe Bose condensation at low momenta l33l l34l . 
However, it is known to neglect important vertex corrections in 
this case @ il as is explained in section II VI 
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Normalized scaling variable: £/B 

FIG. 2. The scaling function /s(£ = t^|p|) = p), 

normalized to the amplitudes A and B, exhibits accurate 
agreement between the nonrelativistic (circles) and relativis¬ 
tic (squares) simulation results. 


notably, the emergence of an effectively conserved par¬ 
ticle number plays a crucial role for the nonequilibrium 
evolution of the relativistic theory. 

The quantitative agreement of the NLO estimates of 
a = 3/2 and /3 = 1/2 for d = 3 with th e full simulation 
results for the relativistic (section III B) and for the non¬ 
relativistic (section |ll B ) theory is remarkable. Though 
we extend the analytic estimates to include a possible 
anomalous scaling in section [VJ we find no strong indica¬ 
tions for a nonvanishing anomalous dimension. Further¬ 
more, we analyze the phenomenon of nonthermal Bose 
condensation from the inverse particle transport towards 
the zero mode. Most importantly, the dynamical genera¬ 
tion of a mass gap m for the (massless) relativistic theory 
from an intriguing interplay of condensation and medium 
effects explains why the relativistic theory exhibits non¬ 
relativistic physics at low momenta |p| < m, as is shown 
in section uni 

The numerically computed fixed point distribution fs 
of equation 0 is given in FIG. [2j Shown is the nor¬ 
malized distribution fs/A as a function of the rescaled 
variable £/B, using the model specific amplitudes A and 
B as defined above. For the relativistic theory, we show 
results for N = 2 field components in FI G. [2] but we 
consider also other values of N in section IIIIl The lat¬ 
tice data for both nonrelativistic and relativistic theories 
turns out to be well described by the fit (dashed line) 


fs(0 


_ A{k> - K<) _ 

(«>-2)(£/B)*< +(2-K < )(£/fl)«> 


(3) 


for £ > 0 with approximately k < — 0.5 and «> ~ 4.5 
for the available range of momenta and times. We note 
that the value for k > is rather close to those obtained 
for stationary turbulence in nonrelativistic systems [19] 
as shown in appendix [A] 


The very good agreement of the exponents and non¬ 
thermal fixed point distributions for the different theories 
in the infrared scaling regime is a striking manifestation 
of universality far from equilibrium. One may use this 
universality, for instance, to learn from experiments with 
ultracold atoms aspects about the dynamics during the 
early stages of our universe [35]. 


II. NONRELATIVISTIC BOSE GAS 
A. Initial conditions 

We first consider a nonrelativistic atomic gas of inter¬ 
acting bosons with s-wave scattering length a(j For a gas 
of density n the average interatomic distance is n -1 / 3 . 
Together with the scattering length a, this can be used 
to define a dimensionless ‘diluteness parameter’ 



For a typical scattering length of, e.g., a ~ 5 nrn and bulk 
density n ~ 10 14 cm -3 the diluteness parameter ( ~ ICC 3 
is very small, and in the following we will always assume 
( < 1. The density and scattering length can also be 
used to define a characteristic ‘coherence length’, whose 
inverse is described by the momentum scale 

Q = VI67ran . (5) 

To observe the dynamics near nonthermal fixed points 
for the interacting Bose gas, an unusually large occu¬ 
pancy of modes at the inverse coherence length scale Q 
has to be prepared. Such nonequilibrium extreme con¬ 
ditions may be obtained, for instance, from a quench 
or nonequilibrium instabilities m More pre¬ 

cisely, for a weakly coupled gas of average density 
n = f d 3 p/(27r) 3 /(|p|) this requires a large characteristic 
mode occupancy 


f(Q) ~ ^ » 1 ■ (6) 

This represents an extreme nonequilibrium distribution 
of modes. The large typical occupancies lead to essen¬ 
tially classical dynamics. Most importantly, the system 
in this overoccupied regime is strongly correlated. 

These properties may be understood from a Gross- 
Pitaevskii equation for a nonrelativistic complex Bose 
field ip: 

id t ip{t,x) = + glV’O,*)! 2 ^ (7) 


2 Here, three spatial dimensions are considered and natural units 
will be employed where the reduced Planck constant (h), the 
speed of light (c) and Boltzmann’s constant (kg) are set to one. 















4 


The coupling g is not dimensionless and determined from 
the mass m and scattering length as g = 4ira/m. The to¬ 
tal number of particles, given by N tota \ = f d 3 x\ip(t, x)| 2 , 
is conserved. 

In the mean-field approximation the effect of the inter¬ 
action term in the Gross-Pitaevskii equation is a constant 
energy shift for each particle, 


A E = 2g(\^\ 2 ) = 2gn = 2gj ^/(|p|), (8) 

which can be absorbed in a redefinition of the chemical 
potential. However, we note that for the very high occu¬ 
pancy © of the typical momentum Q the shift in energy 
is not small compared to the relevant kinetic energy, i.e. 
2 gn ~ Q 2 /2m. Parametrically, this can be directly veri¬ 
fied using ([6]): 

9 I d3 Pf( IpI) ~ 9 Q 3 f(Q) ~ 9%~ ~ 9— 75 ~ — • (9) 
J C mgQ m 

Here we have used that with a = mg/(4n) equation ([5]) 
implies Q = 2 y'mgn and © gives £ = mgQ/(16ir 3l/2 ). 
Most importantly, the energy shift 2 gn is of the order 
of the kinetic energy Q 2 /2m irrespective of the coupling 
strength g. This already hints at a strongly correlated 
system, where the dependence on the details of the un¬ 
derlying model parameters is lost. 


B. Self-similar dynamics from classical-statistical 
simulations 

In order to perform simulations beyond the mean-field 
approximation, we use the fact that the nonequilibrium 
quantum dynamics of the highly occupied system can be 
accurately mapped onto a classical-statistical field theory 
evolution [35l 1391142] . This mapping is valid as long as 
/ 1 for typical momenta, with small enough diluteness 

parameter £ according to I©. In classical-statistical sim¬ 
ulations, one samples over initial conditions and evolves 
each realization according to the classical equation of mo¬ 
tion I©. This equation is solved on a three-dimensional 
grid using a split-step method [23J[24]. Then, observables 
are obtained by averaging over classical trajectories. 

We concentrate on scaling properties of a time- 
dependent occupation number distribution f(t, p). More 
precisely, we consider the two-point correlation function 

f(M',x-x') = + i/)(f',x')f(l,x)) 

( 10 ) 

evaluated at equal times t = t' for spatially homogeneous 
ensembles. Brackets (...) indicate sample averages. In 
spatial Fourier space we define [421 : 

f(t, p) + (27r)V 3) (p)|V’o| 2 (i) = J d 3 xe~ lpx F(t,t,x). 

( 11 ) 



Rescaled momentum: p 


FIG. 3. Rescaled distribution function of the nonrelativistic 
theory as a function of the rescaled momentum for differ¬ 
ent times. The inset shows the original distribution without 
rescaling. 


Because of spatial isotropy, the distribution function de¬ 
pends on the modulus of momentum, and we frequently 
write f(t, |p|). The term ~ \i/j 0 \ 2 (t) coming together with 
the Dirac 4-function represents the condensate part of 
the correlator at zero momentum. The time-dependent 
condensate fraction at zero momentum is then given by 

N o(t) = _ l^o| 2 (*) _ / 12 N 

-^total f d 3 p/(2n) 3 f(t, p) + \il>o\ 2 (t) 


The corresponding expressions for the fin ite volumes em¬ 
ployed will be discussed in section II C| 3 


We consider initial conditions with high occupation 
numbers as motivated in the previous section. Specifi¬ 
cally, we start with a distribution function of the form 


/(0,p) 


vQ(Q- IpI), 


(13) 


which describes overoccupation up to the characteris¬ 
tic momentum Q. The initial condensate fraction is 
taken to be zero, i.e. \ijjo\ 2 (t = 0) = 0, with an ini¬ 
tial /(0,p) = 50/(2 mgQ)Q(Q — |p|). For the plots we 
typically average over 10 realizations on a lattice with 
512 3 points and a lattice spacing a s such that Qa s = 1, 
where we checked insensitivity of our infrared results to 


3 In the quantum theory jlO| denotes the anticommutator expec¬ 
tation value, which has a well-defined equal-time limit, and the 
definition of the distr ibut ion function is obtained by the substi¬ 
tution / —>• / +1/2 in 0- Since our typical occupation numbers 
are large, we drop here the ‘quantum-half’. Furthermore, for the 
class of initial conditions considered, no disconnected part of the 
correlation function arises. This does not exclude a nonzero con¬ 
densate contribution that scales proportional to volume, which 
will be discussed in detail in section III C| 
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cutoff changes. We always plot dimensionless quanti¬ 
ties obtained by the rescalings f(t, p) —> f(t,p)2mgQ, 
t —» tQ 2 /(2m) and p —► p /Q. This reflects the classical- 
statistical nature of the dynamics in the highly occupied 
regime, which has the important consequence that if we 
measure time in units of 2m,/Q 2 and momentum in units 
of Q then the combination f(t, p) 2mgQ does not depend 
on the values of to, g and Q [25] 


The initial mode occupancies (13) get quickly redis¬ 


tributed at the beginning of the nonequilibrium evolu¬ 
tion and then a slower behavior sets in. The latter 
reflects the dynamics near the nonthermal fixed point, 
where universality can be observed. We concentrate on 
the low-momentum part of the distribution and analyze 
its infrared scaling properties. FIG. [3] shows the rescaled 
distribution (t/t re f)~ a f(t, p) of the nonrelativistic the¬ 
ory as a function of (t/t re f)^|p|, where the reference time 
t re iQ 2 /(2m) = 300 is the earliest time shown. The inset 
gives the curves at different times together with the ini¬ 
tial distribution without rescaling for comparison. With 
the appropriate choice of the infrared scaling exponents 
a and /3, all the curves at different times lie remarkably 
well on top of each other after rescaling. This is a strik¬ 
ing manifestation of the self-similar dynamics ([I]) near a 
nonthermal fixed point. The scaling exponents obtained 
are 


a = 1.66 ±0.12, 


/3 = 0.55 ±0.03, 


(14) 


where the error bars are due to statistical averaging and 
fitting errors, which is further described in appendix [Cj 
Comparing these values to |2]), we find that the nu¬ 
merical results (141 agree rather well with the NLO ap¬ 


proximation for a vanishing anomalous dimension, which 
are derived in section El Furthermore, the simulation 
results confirm that a = 3/3 to very good accuracy as ex¬ 
pected from number conservation in the infrared scaling 
regime (see section [Tj|. The values for a and /3 determine 
the rate and direction of the particle number transport, 
since according to ([T]) a given characteristic momentum 
scale K(t 0 ) = K 0 evolves as K(t) = Ko(t/to)~ lJ with 
amplitude f(t, K(t)) ~ t a . Therefore, the positive values 
for the exponents in the infrared scaling regime imply 
that particle number is being transported towards low 
momenta, thus confirming an inverse particle cascade. 


C. Condensate formation 


For the initial conditions (131, there is no condensate 


present at t = 0. However, the inverse particle cascade 
towards the infrared continuously populates the zero¬ 
mode, which leads to the formation of a condensate far 
from equilibrium [24, 2BJ [25] I501 - f52] . Near the nonther¬ 
mal fixed point, the condensation dynamics is expected 
to exhibit scaling behavior and in the following we deter¬ 
mine the relevant scaling exponent. 

Our starting point is the Fourier transformed equal¬ 
time correlation function © with initially vanishing 



FIG. 4. Evolution of the zero-momentum correlation divided 
by volume for the nonrelativistic Bose gas. A power-law ~ t a 
has been fitted in the time interval with self-similar evolution 
for the three largest volumes (dashed lines). Subsequently, the 
results for different volumes converge, signaling the formation 
of a coherent zero mode spreading over the entire volume. 


condensate fraction, |V’o| 2 (^ = 0) = 0. To analyze our 
simulation results for t > 0, we follow references [24J [25] 
and note t hat for a finite volume V the zero-mode contri¬ 
bution in (11) scales with (27r) 3 <5^ 3 ^ ( 0 ) —► V proportional 
to volume. Therefore, we can identify the condensate 
fraction by its scaling behavior as the volume is changed. 
Stated differently, if we divide the correlation function 
by the volume, only correlations which scale with the 
system size are not suppressed at large volumes and the 
condensate fraction is related to the volume-independent 
part. 

FIG.|4]shows the evolution of the zero-momentum cor¬ 
relation V~ 1 F(t, t, p = 0) = F _1 f d 3 xF(t,t,x) for dif¬ 
ferent volumes. These are given in units of Q, ranging 
from VQ 3 = 32 3 to the largest volume ~ 512 3 . Corre¬ 
spondingly, the plotted dimensionless results are rescaled 
as V~ 1 F(t,t,p = 0) —► V~ 1 F(t,t,p = 0)2 mgQ/Q 3 . 
One observes that at early times the evolution depends 
strongly on the volume, as expected in the absence of 
a coherent zero mode spreading over the entire volume. 
However, after a transient regime the zero-momentum 
correlation becomes volume independent. The curves 
corresponding to different volumes converge towards the 
same value, thus signaling the formation of a condensate. 

The double logarithmic plot of FIG. [4] clearly indicates 
that the growth of the zero-momentum correlation pro¬ 
ceeds as a power-law in time. The power-law growth is 
rather well described by the scaling exponent a found 
in (141 from the self-similar evolution of the distribution 
function. These numerical findings are also explained 
in more detail using the analytic scaling solution in sec¬ 
tion m 

Of course, the time needed to fill the entire volume with 
a single condensate increases with volume, which is nicely 
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FIG. 5. Evolution of the condensate fraction for the nonrela- 
tivistic Bose gas for different volumes V. The different curves 
become approximately volume independent after rescaling of 
time by V _1 /“, in agreement with ql5p . 


observed from the data. Using that the parametrically 
slow power-law dynamics dominates the time until con¬ 
densation is completed, we can use the scaling exponent 
a to estimate this condensation time. Taking the value of 
the zero-momentum correlator V~ 1 F{t,t, p = 0) at the 
initial time to of the self-similar regime as U -1 /(to,0) 
and its final value at the time tf as iV’ol 2 (f/), we can 
estimate from V~ 1 F{t,t,p = 0) ~ t a the condensation 
time as 


tf 


to 


/(to,0) 



(15) 


Correspondingly, we define the condensate fraction 
(12) for the case of finite volumes as No/N to ta\ —> 
V^F(t,t,p = 0 )/F(t,t,x = 0), using that N total = 
F(i,t, x = 0). In Fig. [5] we show the evolution of the 
condensate fraction for different volumes as a functi on o f 
the rescaled time t/V 1 '“. Indeed, as predicted by (15), 
the different curves are approximately volume indepen¬ 
dent. One finds that the condensate fraction saturates at 

No/Ntotai ~ 0 . 8 . 


III. RELATIVISTIC SCALAR FIELD THEORY 
A. Initial conditions 

Overoccupied relativistic quantum field theories play 
an important role in early-universe cosmology and high- 
energy collision experiments with heavy nuclei. For a 
large class of inflationary models of early-universe dy¬ 
namics, the accelerated expansion of the universe during 
inflation leads to a large coherent field amplitude of the 
(scalar) inflaton. A subsequent decay via nonequilibrium 
instabilities can lead to highly occupied modes of the 


scalar field with a characteristic momentum Q [53]. A 
similar overoccupation of modes can also occur in heavy- 
ion collisions at ultra-relativistic energies, where a lon¬ 
gitudinally expanding plasma of highly occupied gluon 
fields is expected to form shortly after the collision [551 . 
For expanding systems there are even striking indications 
that the gluon plasma can belong to the same far-from- 
equilibrium universality class as a scalar held theory [T|. 

To be specific, we consider here N real scalar Helds 
ip a interacting via a weak quartic self-coupling A <C 1 in 
three spatial dimensions. The equations of motion for 
the a = 1 ,,N massless Helds are 

- V 2 + g^:<£a(£,x)y> a (t,x)^ <£&(f,x) = 0, (16) 


where a sum over repeated indices is implied. The rele¬ 
vant two-point correlation function we denote as 


F(t, t', X— x') = ^{ip a (t,x)ip a (t' ,x')+<p a (lf 

(17) 

At equal times t = t' this can be used to define a distri¬ 
bution function fit, p) for the relativistic theoryf] 


fit, P) 
u(t, p) 


(27r) 3 A < ' 3) (p)<^>Q (t) = J d 3 xe * px F(f,t,x). 


(18) 

This is in complete analogy to the definition ( [IT] ) for 
the nonrelativistic system and we will refer to the term 
~ (j>o (t) as the condensate part. The only major differ¬ 
ence is the appearance of the dispersion w(t,p) in the 
definition for the relativistic case, which is a consequence 
of the second-order differential equation in time for the 
fields (fl6|). 


Similar to the above discussion for the nonrelativistic 
system, we first characterize overoccupied initial condi¬ 
tions for the typical momentum Q. In a mean-field or 
large-V approximation to leading order, one finds for 
the evolution equation of the correlation function 
in spatial Fourier space (see e.g. [45]): 



d 3 q 

(2tt) 3 


F(t, t, q) 


F{t,t',p) 


0. (19) 


If there is no condensate initially, we can estimate para¬ 
metrically the mean-field correction at sufficiently early 
times as 



, f d 3 p f{t, p) 
J (27t) 3 w(p) 

~ \f{t,Q)Q 2 , 


( 20 ) 


where we have taken a relativistic dispersion w ~ |p| for 
massless particles. One observes that this is of the same 


4 


Similar to the nonrelativistic case, we drop here the ‘quantum- 
half’ as explained in footnote |3l 
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order as the typical kinetic energy term 
the occupancy is as large as 


Q 2 in (191 if 


f(t,Q) ~ ^ > 1 ■ 


( 21 ) 


Common scalar inflaton models for early universe dy¬ 
namics have couplings of order A ~ 10 —13 , such that the 
typical occupancies are extremely large in that context. 
Comparing to ([6]), we note that the dimensionless self¬ 
coupling A plays the role of the diluteness parameter ( in 
the nonrelativistic theory. 


B. Self-similar dynamics from classical-statistical 
simulations 



Rescaled momentum: t* 5 p 


Since the self-similar dynamics can only be observed 
beyond the mean-field approximation, we perform first 
classical-statistical lattice simulations similar to what is 


done for the nonrelativistic theory in section IIB For this 


we solve (16) as classical equations of motion on a three 


dimensional lattice using a leapfrog algorithm mm and 
sample over initial conditions. We extract the occupa¬ 
tion number distribution according to (181 by writing for 

IPl > 0: 


/(C p) =F , t t \ \ZF(t,t',p)d t dt'F(t,t',p) 

W (^P) P \Jd t d t >F(t, t’, p)/F(t, t\ p) t=f 

( 22 ) 

The second equality allows us to identify the dimension¬ 
less distribution (see e.g. [15] ) 


f(t , p) = \/F(t,t / ,p)dtdt>F(t,t , ,p) 


(23) 


The corresponding dispersion relation is then given by 


uj{t,p) = , 


ldtdt'F(t,t',p) 

F{t,t',p) 


(24) 


We emphasize that the notion of particle number is not 
a uniquely defined concept in a relativistic field the¬ 
ory, where total number changing processes are possi¬ 
ble. However, one may always think in terms of the well- 
defined correlation functions appearing in (23). More¬ 


over, this definition turns out to provide an extremely 
useful quasi-particle interpretation even in the strongly 
correlated infrared regime 0, which we will exploit fur¬ 
ther in section (IV). 


Similar to (131 for the nonrelativistic theory, we choose 


overoccupied initial conditions 


1 


/(0,p) ~ j ©(<2- |p|) 


(25) 


with (t = 0) = 0. All quantities shown will be made 
dimensionless by appropriate powers of the scale Q t = 


FIG. 6. Rescaled distribution function of the relativistic two- 
component theory as a function of the rescaled momentum 
for different times. The inset shows the original distribution 
without rescaling. 


yj A e/N, which is obtained from the conserved energy 
density average 

6 = ( \( dt F°)( d tFa) ± ^(dupa^dupa) ± 

' (26) 

where summation over spatial directions i and scalar 


components a is implied. For the initial condition (251, 
one has e ~ NQ A /A such that Q e becomes independent 
of the coupling and the number of field components. 


For the figures we have chosen the amplitude of (25) 
as /(0,p) = 125/A 0(Q — |p|) with Q = 0.8Qe, and 
we always plot rescaled functions f(t, p) —> A f(t , p) and 
F(t, t, p) —I A F(t, t, p) such that these combinations also 
become independent of the coupling. For the N = 2 
component theory, computations were made on a 768 3 
lattice with lattice spacing Q e a s = 0.9 and we averaged 
over five realizations. For N = 4 we employed a 512 3 
lattice with spacing Q e a s = 1.8 and averaged over 18 — 50 
realizations. We checked that for N = 2 all shown results 
are insensitive to the lattice spacing and the volume. For 
the relevant infrared quantities this is to good accuracy 
also the case for the coarser lattices employed for N = 4. 

In FIG. [6] we show the evolution of the occupation 
number distribution for the relativistic N = 2 compo¬ 
nent theory. As in the nonrelativistic case, we plot t~ a f 
against £^|p| to study self-similarity and give the original 
distribution without rescalings in the inset. With appro¬ 
priately chosen exponents, in the infrared all curves lie on 
top of each other after rescaling to remarkable accuracy. 
The measured exponents are 


a = 1.51 ±0.13, 


/3 = 0.51 ±0.04, 


(27) 


and we refer to appendix[C]for details on how we estimate 
the error bars. 
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FIG. 7. Rescaled distribution function of the relativistic four- 
component theory as a function of the rescaled momentum 
for different times. The inset shows the original distribution 
without rescaling. 


FIG. 8. The exponents a and ft with statistical error bars 
for the relativistic N = 4 component field theory extracted at 
different reference times t re f- 


In order to check for a possible dependence of the in¬ 
frared scaling properties on the number of field compo¬ 
nents, we also perform lattice simulations for N = 4. 
FIG. n shows t~ a f as a function t^|p|. The curves corre¬ 
sponding to different times lie again on top of each other 
after the rescalings and we extract the exponents 


C. Generation of a mass gap 


a = 1.65 ± 0.09, P = 0.59 ± 0.03. 


(28) 


These and the universal scaling form of the distribution 
function compare rather well to those for the relativistic 
N = 2 component system as well as the nonrelativistic 
theory. Within the statistical errors we find no indica¬ 
tion for a dependence of the corresponding universality 
class on the number of field components. However, small 
discrepancies in a, P and fg are still possible. These 
could occur in the presence of nonvanishing anomalous 
dimension 77 , which is discussed further in section |V| 

In order to estimate systematic errors, we investigate 
how the values for the exponents a and p depend on the 
reference time t re f at which we start our self-similarity 
analysis. To this end, we perform our analysis for differ¬ 
ent values of t re f and use the method of appendix [C| with 
the distribution function evaluated at several times up to 
t /f re f <4 — 5. In FIG. [ 8 j we show the extracted values 
for a and p as a function of the reference time for N = 4. 
One finds that the mean value of a decreases monotoni- 
cally to about a ss 1.5 while p gets close to a half for the 
transient times at which self-similarity can be accurately 
observed. For the relativistic two-component system, the 
exponents a and P are found to start from comparably 
larger values at early reference times to the ones given 
in (27). We note that the nonrelativistic system of sec¬ 


tion IIB also shows decreasing exponents a and p but 
our runs last not long enough such that their mean val¬ 
ues (141 come as close to the expected values a = 1.5 and 


The values (27) or (28 1 for the exponents of the rela¬ 
tivistic theories agree within errors with the values (fl4|) 


obtained for the nonrelativistic Bose gas. Moreover, they 
are rather close to the analytic results ([ 2 ]) from the large¬ 
ly expansion to NLO for the nonrelativistic theory, while 
they deviate clearly from the corresponding prediction 
assuming a relativistic dispersion as will be shown in sec¬ 
tion [TVl 

Also the universal form of the nonthermal fixed point 
distribution fg accurately agrees for relativistic and non¬ 
relativistic theories in the infrared scaling regime. This 
is shown in FIG. [2] by comparing the two-component rel¬ 
ativistic and the nonrelativistic theory, where fs/A is 
given as a function of t@\p\/B normalized to the nonuni- 
versal amplitudes A and B as described in section |IB| 

To understand the appearance of nonrelativistic dy¬ 
namics, we analyze the dispersion relation w(f, p) accord¬ 
ing to (24), which is shown in FIG. [9] at three different 
times. Although the underlying theory is massless, it can 
be clearly observed that for low momenta the system gen¬ 
erates a mass gap, whereas at large momenta we recover 
a linear dispersion. 

The appearance of an effective mass-like contribution 
can already be understood qualitatively from the approx¬ 
imate evolution equation ( [T9] ) for the correlator modes 
F(t,P, p). In the mean-field approximation, the term 
~ (A/ 6 ) / d 3 q/(2Tr) 3 F(t, t, q) generates a mass-like cor¬ 
rection for the overoccupied initial condition (25). 


p = 0.5 as for the relativistic two-component system. 


To extract the mass gap beyond the mean-field ap¬ 
proximation using the lattice simulations, we fit a time- 
dependent effective mass m(t) from a yjm 2 (t) + p 2 fit to 
the u>(t, p) data. The time evolution of m(t) is shown in 
the inset of FIG. [9] We find that after a quick initial evo¬ 
lution this dispersion relation enters a quasi-stationary 
regime, which is typical for a prethermalized quantity 
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Momentum: p 



FIG. 9. Dispersion relation for the relativistic N = 2 compo¬ 
nent theory at different times for the same parameters as for 
FIG. [6] In the inset, the time evolution of the mass gap at zero 
momentum is shown, which is obtained from a \Jm 2 + p 2 fit 
to <Up at low momenta (+ symbol) and by measuring the oscil¬ 
lation frequency lo c of the unequal-time correlation function 
(x symbol). 


FIG. 10. Evolution of the zero-momentum correlation divided 
by volume for the relativistic N = 4 component theory, simi¬ 
lar to fig. a for the nonrelativistic case. 


D. Condensate formation 


whose transient evolution is governed by an approxi¬ 
mately conserved (particle) number [IB] . 


In addition, we also analyze the oscillation frequency 
oj c of the unequal-time correlation function F(t, 0, p) as a 
function of time t for |p| = 0. Since the zero-momentum 
mode frequency corresponds to the renormalized mass of 
the theory, this provides an independent estimate of the 
mass gap that does not rely on the definition (241 of a 
dispersion. Indeed we find u> c ~ m to very good accuracy 
as shown in the inset of FIG. 1 


In the presence of a mass gap m , low momentum 
modes with |p| < m are expected to behave nonrela- 
tivistically. From FIG. [9] we can estimate the mass to 
be m ~ 0.55Q e ~ 0.69Q for the whole duration of the 
self-similar evolution. We find that this mass scale sepa¬ 
rates rather well the inertial ranges for the inverse parti¬ 
cle cascade towards low momenta from the direct energy 
cascade at higher momenta. This can be observed, for 
instance, from the inset of FIG. |6j where the initial scale 
Q marked by the distribution at t = 0 can be used as a 
reference. 


The emergence of a mass gap in the relativistic the¬ 
ory explains why the dynamics in the infrared regime 
is governed by nonrelativistic physics. Of course, in gen¬ 
eral the presence of a mass gap does not necessarily imply 
universal behavior for sufficiently low momentum modes. 
However, it may be seen as a prerequisite for relativis¬ 
tic theories to belong to the same far-from-equilibrium 
universality class as the Gross-Pitaevskii field theory. 


Since we observe for the relativistic scalar field theory 
the same inverse particle cascade with universal expo¬ 
nents and scaling function as for the nonrelativistic sys¬ 
tem, one expects condensation to proceed with the same 
exponent a as found in section IIC However, this is not 
entirely trivial since scatterings off the condensate play 
an important role in the inertial range for the direct en¬ 
ergy cascade towards higher momenta, where both theo¬ 
ries do not belong to the same universality class [ 3111 - 


In order to clarify this, we analyze the growth of the 
condensate during the self-similar regime for the N = 4 
component theory along the same lines as in section [TTC] 
The results are given in FIG. [To] for different volumes, 
ranging from V= 58 3 to the largest volume ~ 922 3 . 
Shown is the time evolution of the zero-mode of the cor¬ 
relation function (181 divided by volume. We compare 
the curves to the expected power-law behavior ~ t a with 
the exponent a obtained before from the self-similarity 
analysis. As in the nonrelativistic case, the dynamics is 
well described by such a power-law until the results be¬ 
come volume independent, thus signaling the formation 
of a coherent condensate 4% f° r the entire volume. 


We emphasize that the observed power-law is restricted 
to the formation of the condensate during the transient 
self-similar regime, where we find the universal exponent 
a to govern the dynamics. In particular, we do not dis¬ 
cuss here the subsequent late-time approach to thermal 
equilibrium, where total particle number changing pro¬ 
cesses in the relativistic theory can make important dif¬ 
ferences as compared to the number conserving nonrela¬ 
tivistic system [21] . 
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IV. NONTHERMAL FIXED POINTS FROM 
VERTEX-RESUMMED KINETIC THEORY 

A. Self-similarity 


with a generic ‘collision integral’ C[f](t, p), which de¬ 
pends on the theory and the approximation employed. 
For the self-similar distribution (29), the scaling behav¬ 
ior of the collision integral is then given by 


While many aspects of transport of conserved charges 
associated to wave turbulence in kinetic theory have long 
reached textbook level, there is still rather little known 
about the analytic description of turbulent behavior in 
nonperturbative regimes of quantum field theories. 

Perturbative kinetic theory has been successfully ap¬ 
plied to the phenomenon of weak wave turbulence [3, 
[2ni mi- It has also been employed in the literature 
to describe infrared phenomena such as Bose condensa¬ 
tion [33j, 33] ■ However, in this case the approach ne¬ 
glects important vertex corrections since the large occu¬ 
pancies at low momenta lead to strongly nonlinear dy¬ 
namics Here we apply a vertex-resummed kinetic 

theory, which is based on an expansion in the number of 
field components N to next-to-leading order psissiezi. 
It extends well-known kinetic descriptions [31 HOI HI] to 
the nonperturbative regime of overoccupied modes. This 
allows us to gain analytic understanding of the formation 
of a dual cascade, in which particles are also transferred 
towards low momenta leading to Bose condensation. 

Similar nonperturbative descriptions have been em¬ 
ployed before for the scaling behavior of stationary trans¬ 
port, which is time-translation invariant [2 US El]. Since 
our isolated systems out of equilibrium cannot realize 
stationary transport solutions in the absence of external 
driving forces, we have to consider here the more gen¬ 
eral notion of a self-similar evolution. Though this is 
not time-translation invariant, the dynamics in this case 
is described by time-independent scaling functions and 
scaling exponents. 

A self-similar evolution of the distribution function 
f(t, p) for a spatially homogeneous and isotropic system 
is characterized as 


C[f](t,p) = s M C[/](s 1/l3 t,sp) = t P ^C[f s ]{ M^p), 

( 31 ) 

where /i is a function of scaling exponents such as a and 
0. Substituting this scaling form into the kinetic equa¬ 
tion leads to the time-independent fixed point equation 
for fs( P): 


[a + 0p- V p ]/s(p) = C[/ S ](l,p), (32) 

and the scaling relation: 

a-l = -dp ,. (33) 

This follows from comparing the LHS of the kinetic equa¬ 
tion, 

[t a fs(t p p)] = t a ~ 1 [a + 0q- Vq]/s(q)| q=t/3p , (34) 


to its RHS given by ([31]) . 

Further relations can be obtained by either imposing 
energy conservation or particle number conservation if 
applicable. For constant 


n = 


(35) 


one obtains the relation for 

particle conservation: a = 0d. 
Similarly, with 


(36) 


f(t,p) = s a /f > f(8-Wt,sp) (29) 

with the real scaling parameter s and exponents a and 
0. Again, all quantities are considered to be dimension¬ 
less by use of some suitable momentum scale. Choosing 
s-Wt = 1 we recover (Jll) , where the time-independent 
scaling function fs(t^ p) = f(l,t^p) denotes the fixed 
point distribution. This scaling form represents an enor¬ 
mous reduction of the possible dependence of the dynam¬ 
ics on variations in time and momenta, since t~ a f(t, p) 
depends on the product |p| instead of separately de¬ 
pending on time and momenta. Therefore, an essential 
part of the time evolution is encoded in the momentum 
dependence of the fixed point distribution /s(p). 

For the time evolution of the distribution function 
f(t, p) we write 

[/]((, p) (30) 


t = j Wf “ (p)/(f ' p) = / Wf “ <q)/s(q> 

(37) 

one obtains from 

energy conservation: a = 0{d + z ), (38) 

where the dispersion is taken to scale with the dynamic 
exponent z as 

w(p) = s -2 w(sp) (39) 

with z = 1 for the gapless relativistic and z = 2 for the 
nonrelativistic theory. One observes that there is no sin¬ 
gle scaling solution conserving both energy and particle 
number. As outlined already above, in this case a dual 
cascade is expected to emerge such that in a given iner¬ 
tial range of momenta only one conservation law governs 
the scaling behavior. 
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B. Perturbative scaling behavior 


We first review some perturbative results for later com¬ 
parison. In perturbative kinetic theory, when two parti¬ 
cles scatter into two particles, the collision integral is of 
the form 

C 2 " 2 [f](t,p) = J dU 2 ^ 2 (p,l,q,r) 

X [(/p + 1)(/1 + l)/ q /r - /p/i(/ q + l)(/r + 1)] , (40) 


where we write /(Z, p) = / p suppressing the global time 
dependence to ease the notation. The details of the 
model enter f dU 2o2 (p, 1, q, r), which for the example 
of the relativistic iV-component scalar field theory with 
quartic A/(4!ZV)-interaction of section III reads: 


<m 2 


’(p,l,q,r) = A 


2 N + 2 f d d l d d q d d r 


6 N 2 J ( 2n) d ( 2n) d ( 2Tr) d 

x (27T) d+1 6^ (p + 1 - q - r) 
v ; VH-r M ; 2w p 2wi2o; q 2w r 

(41) 


with w p = y/p 2 + m 2 . 

The expression represents a standard Boltzmann equa¬ 
tion for a gas of relativistic particles, which is not ex¬ 
pected to be valid if the occupation numbers per mode 
become too large. Parametrically, for a weak coupling A 
a necessary condition for its validity is f p <C 1/A since 
otherwise strongly nonlinear effects become significant as 
will be explained in detail in section [IV C| On the other 
hand, scaling is expected for not too small occupation 
numbers per mode, which we discuss now. For the cor¬ 
responding regime l€/ p Cl/A one may use the above 
Boltzmann equation, which approximately becomes 


Apart from the 4-vertex interaction considered, it will 
be relevant to investigate also scattering in the presence 
of a coherent field such that an effective 3-vertex appears. 
To keep the discussion more general, we may write for the 
scaling behavior of a generic collision term for Z-vertex 
scattering processes 

C {1) [f] (t, p) = C® [/] sp), (45) 

where the scaling exponent 

m = (Z — 2 )d — (Z + 1) — (Z — l)a//3 (46) 


follows from similar arguments as exemplified for the 
4-vertex interaction. 

Using the scaling relation (33) and particle conserva¬ 
tion (36) gives the perturbative solution for 

d 


rel. particle transport: a = — 


l + 1 


•e—m- (47) 


Similarly, with (|38|) for relativistic z = 1 one finds the 

d- 


perturbative solution for 
rel. energy transport: a = — 


1 


1 


21 - 1 ’ 


21-1 


■ (48) 


For instance, for quartic self-interactions the perturba¬ 
tive energy transport is characterized by a = — (d + l)/7, 
and /3 = —1/7, where the latter is independent of the 
dimensionality of space d. Likewise, for an effective 
3-vertex in the presence of a coherent field one has for 
the energy transport a — — (d + l)/5 and /3 = —1/5. In¬ 
deed, the latter values for the scaling exponents describe 
well the energy transport at higher momenta of the 
relativistic scalar field theory for d = 3 of section 00 - 
In particular, their negative values indicate the direction 
of the transport from lower to higher momenta. 


t-/p - j dU 2o2 (p,l,q,r) 

X [(/p + /0/ q /r - /p/l(/ q + fr)} . (42) 

Apart from the energy density e also the total particle 
number density n is conserved by the collision term. 


For the relativistic theory, the scaling assumption (29) 
should be valid for sufficiently high momenta |p| m 
such that the dispersion is approximately linear with 
w p ~ |p|- In this case one obtains for the scaling of the 


collision integral (41) and (42) of the theory with quartic 
self-interaction: 

C 2 ~ 2 [f](t, p) = sp) , (43) 

where the scaling is described by 

Pi = (3d — 4) — (d + 1) — 3a/(3 = 2d-5-3a//3. (44) 

The first term in brackets comes from the scaling of 
the integral measure, the second from energy-momentum 
conservation for two-to-two scattering and the third from 
the three factors of the distribution function appearing 
in p2t). 


We now turn to the nonrelativistic theory of section [TT] 
The perturbative kinetic equation for 2 <->■ 2 scattering is 
given again by (42), however, with the nonrelativistic 


dO n r (pA,q,r) ~ g 


d d l d d q 


d d r 


(2ir) d (2jr) d (2ir) d 

x (27r) d+1 (p + 1 — q — r) 5(w p + wi — w q — u r ) (49) 

in the absence of a condensate for w p = p 2 /(2?n). The 
scaling analysis follows along the same lines as before, 
but without the inverse-frequency factors from the rel- 
ativistically invariant measure appearing in (41). Ac¬ 


cordingly, generalizing again to Z-vertex scatterings, one 
obtains for the scaling relation (331 in the nonrelativistic 
case: (l — 2)a =/3[(l — 2)d—2] — 1 . This leads with (36) 
and (38) to the solutions for 


nonrel. particle transport: a = (3 = — (50) 

and 

, d+2 n 1 

nonrel. energy transport: a = -jy, p = -^y- 

( 51 ) 
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For instance, for quartic self-interactions the perturbative 
particle transport would be described by a = —d/2, and 

P = — 1 / 2 . 

We emphasize that all the above perturbative esti¬ 
mates with negative values for a and /3 cannot account 
for the inverse particle transport observed from the sim¬ 
ulations for d = 3 in sections IIB and IIIB Of course, 
perturbation theory is not expected to be applicable to 
the overoccupied infrared modes and one has to employ 
an alternative description, which we consider next. 


C. Scaling behavior in the overoccupied regime 


Remarkably, the overoccupied regime can still be de¬ 
scribed in terms of a generalized kinetic theory by taking 
into account vertex corrections. For the iV-component 
field theory, these corrections can be systematically com¬ 
puted from an expansion in the number of field compo¬ 
nents N to next-to-leading order (BJ [35) 137] . The NLO 
corrections take scattering events up to infinite order into 
account. This allows us to describe even strongly corre¬ 
lated regimes, where the typical mode occupancies © or 
(21) are inversely proportional to the diluteness or cou¬ 


pling parameter, respectively. 

Effectively, the change to the perturbative kinetic 
equations (40) or (42) is the appearance of a time and 


momentum dependent matrix element squared: 


A 




(52) 


More precisely, one finds at NLO of the expansion in the 
number of field components for the relativistic theory a 
kinetic equation where (41) is replaced by: 


f mt nr ci/ , x 1 f d d l d d q d d r 

JiH [/](<, P,l,q,r)=—y^^ Fwi 

x p*)" ■ «<■%+1 - q - r) ^rrrr"' 1 

ZLUpZOJiZLOqZLUr 


X Aeff[/](t,p,l,q, r), 
with dispersion w p = \/p 2 + m 2 . Here 


(53) 


A 2 

Aeff(*,P,l,q,r) = — 


1 


|i + n fl (t,w p + wi,p + i )| 2 

l 

+ |1 + n fl (f,Wp — Wq, p — q )| 2 
+ |l-f-n- R (t,o;p-w r ,p-r)| 2 _ *' 54 ' ) 

incorporates the vertex corrections for the different scat¬ 


tering channels. They are depicted in FIG. 11 and may 
be viewed as coming from an effective interaction, which 
involves the exchange of an intermediate particle whose 
four-momentum equals p+l , p — q and p—r, respectively. 



FIG. 11. Illustration of different scattering channels. The 
vertex correction at NLO may be viewed as an effective inter¬ 
action, which involves the exchange of an intermediate parti¬ 
cle. Left: The incoming particles with momenta p and l join 
into an intermediate particle that eventually splits. Middle 
and right: The incoming particle with momentum p emits 
the intermediate particle and becomes the final particle with 
momentum q or r, respectively. 


The appearance of the ‘one-loop’ retarded self-energy, 


n l (f,w, p) = — 


A f d d q f(t, p —q) 


12 J ( 2 ir) d WqWp-q 

1 1 


H - ^p—q ^ ^q q ^ 

1 1 


CJq — CJp_q + UJ CJq + CJ p _q + UJ 


(55) 


in the denominator is the result of a geometric series 
summation of the infinite number of NLO processes. It 
should be emphasized that n K , and thus also A 2 ff , is time- 
dependent because it is expressed in terms of the evolving 
distribution function. The above expressions correspond 
to the on-shell limit of the evolution equations presented 
in reference [B: and their relation to the underlying field 
theory is further discussed in appendix |B 3[ 


For sufficiently high momenta, the self-energy (55) be¬ 
comes small such that n fl (i, w,p) <C 1 in the squared 
matrix element (54) and we recover the perturbative ex¬ 
pression A 2 ff 


-¥ X 


In this case, the kinetic equation 
corresponds to ( |41| with the prefactor for large N. Since 
the scaling solutions of section [IVB] do not depend on N , 
one gets the same results in the high momentum regime. 

For low momenta the case n- R (t,w,p) 1 can be¬ 
come relevant, which changes the situation dramatically. 
This becomes even more involved if a mass gap exists, 
as seen in the simulations of section IIII Cl An effective 
mass gap is typically expected because of medium effects 
even if the mass parameter of the underlying microscopic 
theory is set to zero. In that case the infrared modes be¬ 
have effectively nonrelativistically, which allows one to 
observe the same scaling behavior of the relativistic and 
the nonrelativistic theory in this regime. 

The corresponding kinetic equation for the nonrela¬ 
tivistic theory may be obtained from the relativistic case 
(54), assuming that the mass appearing in the disper¬ 


sion relation is much larger than the typical momenta. 

= y/i 


Expanding uj p = i/p 2 + m 2 ~ m + p"/2 m and insert¬ 
ing this into (53), we have <5(w p + uj\ — w q — w r ) —► 
<5((p 2 + l 2 — q -r 2 )/(2m)) and 2w p 2wi2o;q2w r — > 16to 4 
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to lowest nonvanishing order. Then g ~ A/m 2 leads to 
the nonrelativistic form of the equation. Alternatively, 
one can obtain the corresponding result from a nonrel¬ 
ativistic TV-component complex scalar field theory jT5j 
by performing again the 1/iV expansion to NLO. This is 
further described in section El 

Accordingly, we consider for the nonrelativistic case a 
kinetic equation with 


dttn^°[f](t, P,l,q, r) 


d d l d d q d d r 


(2n) d (2ir) d ( 2Tr) d 
x (27r) d+1 4 (d) (p + 1 - q - r) S(uj p + ui\ — w q — w r ) 

x 5eff [/](*. P,q)» (56) 

where uj p = p 2 /(2m) and 


5eff(^P>q) = 


H-HS^Wp-Wq.p-q) 


(57) 


incorporates the vertex corrections. The compact form 
of (57) is achieved by using the symmetries of (56) and 


the nonrelativistic ‘one-loop’ retarded self-energy reads 


n 


(f,o;,p) = 2 g J p-q) 


Wq — W D _a — w 


Wq - Wp-q + 


(58) 


By comparison to its relativistic counterpart, this may 
also be obtained directly from (55) if evaluated as in 


(571 by expanding the relativistic dispersion for small mo¬ 


menta and taking the dominant contributions, where the 
constant mass term cancels in the respective frequency 
sums. 


From the self-similar behavior of the distribution (291, 
we can deduce the scaling property 


nnr(^ p ,p) = s a/ld d+2 n*(s 1//3 t,uj sp ,sp). (59) 


Since for the relevant cases (36) and (38) we have 
a/P > d, we conclude by keeping sp fixed that n„ r can 
become large in the infrared. In this case, we can use 
n^ r (t,o;p, p) 1 to find the scaling behavior of 

9 %s(t, p, q, r) = s -2(«//3-d+2) S p, sq , sr ). 

(60) 

Therefore, we expect the effective matrix element squared 
to become small in the infrared in accordance with re¬ 
lated studies |2j. In turn, we will see in the following 
that at the same time the distribution function f(t, p) 
can grow significantly, which results in a ‘Bose enhance¬ 
ment’ of scatterings that counteracts the diminished ef¬ 
fective coupling. 


We use (561 in the corresponding kinetic equation (42) 


and find for the scaling of the collision term 

C L °[/](i, P) = 

= t a -^C^°[f s ]( M^p 


Remarkably, with this the exponent P can be obtained 
solely from the scaling relation (33) without using in ad¬ 


dition energy or particle conservation, whereas the dif¬ 
ferent solutions for a arise from imposing (36) or (381, 
respectively. We thus find in the overoccupied infrared 
regime for 


1 ( particles: a = d/2 

nonrel tmn ’ p0Tk 11 = 2 ° f \ energy: a = (d + 2)/2 

(62) 

This is a central analytic result of this work. In contrast 
to the previously known negative scaling exponents 
from perturbative estimates given in section |IVB[ one 
observes that the positive values of a and P obtained 
here describe an inverse particle transport with growing 
occupation number in the infrared. The quantitative 
agreement of the NLO estimates a = 3/2 and /3 = 1/2 
for d = 3 with the full simulation results of section IIII Bl 
for the relativistic and section Hi Bl for the nonrelativistic 
theory is remarkable. Both the approximate analytic 
and the full simulation results within their numerical 
accuracy indicate no strong dependence on N. 

For comparison, we finally also analyze the relativistic 
kinetic equation with (531 in the absence of any mass 


gap. To this end, we use u p = |p| and proceed as in the 
nonrelativistic case to find scaling relations for a and p. 
With 

C NLO [f}(t lP ) = t a ~P C' NLO [ /s](l,t /3 p) (63) 

we find for the 

f particles: a = d 

rel. transport: P = 1 of < , (64) 

( energy: a = d+ 1 v ' 

These estimates indicate that the simulation results of 


section IIIB for the relativistic theory cannot be inter¬ 
preted in terms of massless scaling in the infrared, which 
is explained by the presence of a mass gap in section[Til C| 


V. ANOMALOUS SCALING 

The kinetic description of section [IV| assumes canonical 
values for the dynamic exponent z. which describes the 


scaling of the dispersion (39). For the relativistic theory 
without a mass gap, u: p =/p| is employed. For the Gross- 
Pitaevskii theory, in the presence of a condensate, the 
approximate (Bogoliubov) dispersion is given by [47, 




(65) 


At larger momenta, or in the absence of a condensate, 
one recovers u> p = p 2 /(2m), while for low momenta one 


(61) has w p = c|p| with c = ^g\ipo\ 2 /m. 
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From a field theoretic point of view, these dispersions 
with integer-valued z are implemented using a canonical 
spectral function, such as the free-held form 

Po(p° i P) = 27rsgn(p°) S (( p °) 2 - w 2 ) (66) 

for the relativistic theory. In principle, nonperturbative 
scaling phenomena may involve an anomalous scaling ex¬ 
ponent for the full spectral function p(p°, p) of the inter¬ 
acting theory. Using spatial isotropy we may write 

p(p°, P) = s 2 ~ v p{s z p°,sp) , (67) 

with a nonequilibrium ‘anomalous dimension’ 77 . The dy¬ 
namic scaling exponent 2 appears since only spatial mo¬ 
menta are related by rotational symmetry and frequen¬ 
cies may scale differently also in the relativistic theory 
because of medium effects. 

Small discrepancies between the results of our full nu¬ 
merical simulations and the analytic estimates of previ¬ 
ous sections could possibly be rooted in the canonical 
assumption of an integer-valued z and rj = 0. This con¬ 
cerns mainly the infrared regime, where strongly nonlin¬ 
ear behavior occurs. For instance, corresponding infrared 
scaling phenomena in thermal equilibrium near continu¬ 
ous phase transitions can exhibit non-trivial scaling expo¬ 
nents with a noncanonical 2 and a nonzero (but typically 
small) value for 77 , which is also captured by the NLO 
approximation employed [48]. 

Therefore, we consider in this section a field theoreti¬ 
cal calculation of the self-similar behavior near nonther- 
mal fixed points taking into account the possibility of 
anomalous scaling. It is again based on the two-particle 
irreducible (2PI) generating functional in quantum field 
theory, which is expanded up to next-to-leading order in 
the number of field components N OH HZ]. However, 
without using the additional assumption of a canonical 
form for the spectral function underlying the kinetic the¬ 
ory of section [TV] 


A. Nonrelativistic field theory 


in order to have a compact notation for the four different 
two-point functions that can be built from two complex 
fields. 

The other set of linearly independent two-point func¬ 
tions can be conveniently expressed in terms of the ex¬ 
pectation value of the anticommutator of two fields as 

F a b(x,y) = l l ({i> a (x), , 4>t(y)})c, (71) 


where the subscript ‘c’ refers to the connected correla¬ 
tor. While in thermal equilibrium the anticommutator 
and commutator expectation values are related by the 
fluctuation-dissipation relation, in general they are lin¬ 
early independent for systems out of equilibrium 0 The 
absence of a fluctuation-dissipation relation is a crucial 
property of the scaling behavior near the nonthermal 
fixed points discussed. 

The evolution equations for F ab (x,y) and p a b(x,y) at 
NLO in the 2PI 1/N expansion are known [TSJ [Ml [57] 
and given in appendix [B] for completeness, where we 
also outline the employed gradient expansion to lowest 
order. It is important to notice that, because of the 
lowest order gradient expansion, the spectral function is 
time-independent according to (B7l. For the spatially 
homogeneous system, it is convenient to Fourier trans¬ 
form the two-point functions with respect to their rela¬ 
tive coordinates x — y and to define the ‘time’ variable 
as t = {x° +y°)/2. One obtains, using ( F p ) ab = F ab (t,p) 
and ( p p )ab = Pab{p ) as a compact matrix notation in 
(a, 6 )-index space: 


x <5 {d+1) ip - q + l - r) g 2 s [F\(t,p- q) 
x | 2 TV [<y s F p F q ] Tr [Fj p r ] 

+ Tr [u 3 F p pq\ Tr [Fj F r ] 

— Tr [a 3 p p F q ] Tr [F) F r ] j, (72) 


For the nonrelativistic scalar field theory we consider 
a single complex scalar field ip(x) with the action 

S[ip, t/>*] = J d d+1 x <|V (id x 0 + V - | W *) 2 j • 

( 68 ) 

In the classical approximation this gives rise to the Gross- 
Pitaevskii equation ([T]). In the corresponding quantum 
theory, the spectral function is given by the expectation 
value of the commutator of two Heisenberg field operators 
as 

Pab{x, y) = ([ tp a (x ), V'l(y)]) • (69) 


where Tr[F p ] = F aa (t,p) and a 3 = diag(l, —1) denotes 
the third Pauli matrix. The time- and momentum- 
dependent effective coupling squared corresponding to 
(571 reads 


9es(t,p) = 


9 


|i + n«(t,p )| 2 

with the retarded self-energy 
d d+1 - 


(73) 


/ d d + 1 n 

(2tt )d+i Fab ( t ’ q -P) G ba(9 ) (74) 


Here the index notation for a,b = 1,2 employs 


(70) 


5 For an introductory text, see e.g. reference m- 
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in terms of the retarded propagator G(/ b (p) as is further 


discussed in appendix B 3 


To make contact with the definition of the distribution 
function in ( |11[ ), we note that with the notation F = 
F aa /2 this can be written for p ^ 0 as: 

/(*) P) + \ = \ J i^F aa (t,p 0 , p), (75) 

where we write the ‘quantum-half’ for completeness 
though we always consider high typical occupancies such 
that it can be neglected. With the field theoretical def¬ 
inition of the distribution function at hand, we can ob¬ 
tain a corresponding ‘collision term’ for d.fjt , p )/dt = 
C[F](t,p) from the p°-integral of the RHS of (72). 

Applying (29) to the distribution function with 


fit, P) 1 for typical momenta, self-similar behavior 
for the correlators is described by 

F ab (t,p°, p) = s Q /^F a6 ( S -^MV, S p) (76) 


in addition to the scaling behavior (67) for the spec¬ 
tral function p a b(p°, p)- This means Fab(t,P°, P) = 
t a+pz F Si ab{t Zp P°,t^p) with F S ,ab(p°, P) = F ab {l,p°,p). 
In particular, f s { p) = / dp 0 /(2-K)F s ^ aa (p°, p)/2 in this 
highly occupied scaling regime according to ([75]) . With 


g 2 eS [F](t,p°, p) = t-^+Pe-v-d)] g ^[F s ](t^p°y p) 

(77) 

one finds 


C[F](t, p) =t a ~^ 2 ~^ C[F s \{tP p). 


(78) 


We can now proceed in complete analogy to the anal¬ 
ysis of section |TV| from which we find the solution to the 
scaling relation: 





(79) 


where the dimensionality d and the exponent a have 
dropped out. In addition we have the time-independent 
equation for the nonthermal fixed point function: 


\a + /3p- V p ]/s(p) =C[F s ](p) 


(80) 


corresponding to (32). 


To obtain a second scaling relation for the determina¬ 
tion of a, it is important to note that the particle number 
density n = f d d+l p/{2n) d+1 F aa {t,p)/2 and the energy 
density e = f d d+1 p/(2n) d+1 p° F aa (t,p)/2 are conserved 
during the evolution. We find: 


d 

(particle transport), 

(81) 

2-77 
d + z 
~2 — 77 

(energy transport). 

(82) 


Of course, taking rj = 0 and z = 2 we recover the re¬ 
sults of section |IV C| Remarkably, one observes that the 


scaling exponent for the nonrelativistic particle cascade 
is independent of the dynamic exponent 2 . This has the 
important special consequence that the same scaling be¬ 
havior for particle transport is found no matter whether 
the low-momentum dispersion is quadratic, or linear as 
for the Bogoliubov dispersion (651 in the presence of a 
condensate. 

Comparing the expressions for particle transport to our 
simulation results of section (IIB), within the numerical 


errors we see no strong indication for any deviation from 
the canonical value 77 = 0 employed for the analytical 
estimates of section El 

We may also consider the growth of the condensate 
during the self-similar evolution, which we investigated 
numerically in sections [TT| and III With F ab (t,t,p) ~ 
f dp 0 F ab (t,p° , p) and using ( 76) we obtain the scaling of 
a characteristic mode of the equal-time correlator ~ t a . 
Since a > 0 for not too large 77 , the condensate zero¬ 
mode is expected to grow as a power-law in time. The 
value a ~ 3/2 for d = 3 rather accurately describes our 
numerical findings. 

For completeness, we indicate how to recover the ki¬ 
netic equation of section |IV C| from the evolution equa¬ 
tion (72) for the anticommutator expectation value. 


First, one inserts the free-held form of the spectral func¬ 
tion uni 


P(p) = 2tt 


(S(p° — 1 

{ 0 


; P ) 0 

—S(p° -t 


J P ) 


(83) 


where w p = p 2 /2m. Second, one can define an off-shell 
distribution function f(t,p) as 


F a b{t,p) = ( f(t,p) + - ) p ab {p) 


(84) 


with — f(t,—p ) = f(t,p) + 1 in accordance to (75). In¬ 
serting all this yields for f(t, p) 1 the nonrelativistic 
kinetic equation employed in section pV C| 


B. Relativistic field theory 


We now consider the relativistic quantum field the¬ 
ory for a Af-component real scalar field p a (x) with a £ 
{1,..., A^} in d spatial dimensions. Its classical action is 
given by 


S' = / d d+1 


: | ^VadptPa - ^ (^a^a) 2 j • 

(85) 


The corresponding classical equation of motion for m = 0 
is the Klein-Gordon equation given in (16). 

Following the discussion of section |IV C| we will focus 
here on the case without a mass gap, since otherwise 
the above nonrelativistic results apply. Along the lines 
of the previous section, we introduce spectral and sta¬ 
tistical functions as the commutator and connected an¬ 
ticommutator expectation values of two Heisenberg field 
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operators. In the following we will assume that the corre¬ 
sponding F ab = F Sab and p ab = p 5 ab . According to ( [l8| ) 
and taking into account the quantum-half, we get the 
distribution function for sufficiently high occupancies for 
p ^ 0 as: 


f(t, p) 


dp 0 
27r 


2 P°F(t,p). 


( 86 ) 


The derivation of the evolution equation for /(f, p) fol¬ 
lows along the same lines as for the nonrelativistic theory 
of the previous section. Using the 2PI 1/N expansion to 
NLO, the lowest-order gradient expansion leads to 


df 1 f c 

di (t,P) ~ 6iV(27r) 2d + 3 J 0 


dp 0 dl°dq°dr 0 I d a ld a qd a r 


di jd. 


x 6 (d+1 \p + l — q — r) A l s [F](t,p, l,q,r ) 

(p P Fi + F p pi)F q F r - FpFi ( p q F r + F q p r ) , (87) 

where F p = F(t,p), and A 2 ff is defined as 
A 2 

= y [v e s(t,P + l) 

+ v e ff (t,p-q) + v e s(t,p-r)] (88) 

with the vertex function v e g(t,p ) given in ( B22[ ). 

With the above kinetic equation we now search for self¬ 
similar solutions (291. For this we have to write 


p) = s a /P+ 2z Fis-^t, s z p°, sp), 


(89) 


in addition to the scaling behavior (671 for the spectral 
function. Inserting the scaling behavior into the evolu¬ 
tion equation ((871) leads to the solution 


1 


2 — rj — z 


(90) 


and an equation for the universal scaling function corre¬ 
sponding to (80). Imposing particle number and energy 
conservation yields 

d 


a = 


a = 


2 — rj — z 
d + z 


(particle transport), (91) 

(energy transport). (92) 


2 — ii — z 

For p = 0 and z = 1 these values agree with the results 
in the absence of a mass gap of section IV C We finally 
note that the relativistic kinetic equation of that section 
can be obtained from (87) by using the corresponding 
definition (84) with the free spectral function (66). 


VI. CONCLUSION 

A consistent picture has emerged for the universal self- 
similar dynamics of relativistic and nonrelativistic field 
theories near nonthermal fixed points. The results of full 


numerical simulations are well described by analytic re¬ 
sults based on the vertex-resummed kinetic theory. The 
latter extends well-established kinetic descriptions to the 
nonperturbative regime of overoccupied modes. 

The vertex-resummed kinetic theory links the pertur¬ 
bative phenomenon of weak wave turbulence relevant at 
higher momenta to the nonperturbative physics of strong 
turbulence and vorticity of the underlying field config¬ 
urations in the highly nonlinear infrared regime. It is 
striking that the range of validity of kinetic descriptions 
can indeed be extended to capture these very different 
regimes. 

For the examples of nonrelativistic (Gross-Pitaevskii) 
and relativistic scalar field theory with quartic self¬ 
interactions, we have seen that the infrared scaling ex¬ 
ponents as well as scaling functions agree. This becomes 
possible because of the emergence of a mass gap in the 
relativistic theory. In contrast to the previously known 
negative values for the scaling exponents a and /? from 
perturbative estimates, we find that their positive val¬ 
ues a = d/(2 — rj) and fi = 1/(2 — rj) obtained for small 
anomalous dimension ry describe an inverse particle trans¬ 
port with growing occupation number in the infrared. 
The growth exponent a is found to describe also the far- 
from-equilibrium formation of the Bose condensate. 

Moreover, this nonrelativistic particle transport solu¬ 
tion has the remarkable property to be independent of the 
dynamic scaling exponent z. As a consequence, this so¬ 
lution applies equally well to a dispersion with quadratic 
momentum dependence or a possible linear behavior be¬ 
low the characteristic coherence momentum scale in the 
presence of a Bose gas condensate. 

The corresponding dynamic universality class turns 
out to be remarkably large, encompassing both relativis¬ 
tic as well as nonrelativistic quantum and classical sys¬ 
tems. As a consequence, the applications can range from 
table-top experiments with ultracold quantum gases to 
inflationary dynamics during the very early stages of our 
universe. 
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Appendix A: Stationary transport 

For comparison, we summarize in this section the anal¬ 
ysis for stationary transport. This is a time-translation 
invariant problem, which can be realized in the presence 
of suitable sources and sinks. In contrast, the isolated 
systems out of equilibrium considered in the main part 
of this work are characterized by a self-similar evolution, 
which is not time-translation invariant. 

Stationary transport in perturbative kinetic theory has 
been intensively studied in the context of weak wave 
turbulence [3 m IS]. Similarly, the nonperturbative 
regime of strong turbulence has been investigated in de¬ 
tail and the following discussion is based on the litera¬ 
ture SI HUGS, where we point out the relevant aspects 
for the comparison to the case of self-similar evolutions. 


Stationary wave turbulence is characterized by a scale- 
independent flux A(k), for which the respective integral 
does not depend on the integration limit k. To this end, 
one considers scaling solutions 

fp = & fsp ; k-’p = $ ^sp 5 (A4) 


with occupation number exponent k and assuming a lin¬ 
ear dispersion relation relevant for momenta |p| 3> m. 
Since the physics is scale invariant, we can choose s = 
l/|p| such that / p = |p| -K fi and u p = |p| uq. 

Using these scaling properties, one obtains for the col¬ 
lision integral (41) and (42) of the theory with quartic 
self-interaction: 


C 2 o 2 (p) = s~^C 2 ^ 2 (sp), 
where the scaling exponent is given by 


(A5) 


1. Perturbative regime: Weak wave turbulence 


= (3d - 4) — (d + 1) — 3k = 2d - 5 - 3k . (A 6 ) 


In perturbative kinetic theory, when two particles scat¬ 
ter into two particles, the time evolution of the distribu¬ 
tion function f p (t ) for a spatially homogeneous system 
is given for not too small occupation numbers by (421. 


For the example of the relativistic iV-component scalar 
field theory with quartic A/(4!IV)-interaction this is fur¬ 
ther specified by (41). The phenomenon of weak wave 


turbulence is expected for occupation numbers per mode 
in the regime 1 -C f p <C 1/A, where that perturbative 
description is expected to be valid. In section |A~2| below, 
we consider the nonperturbative case of an overoccupied 
system. 

For number conserving 2o2 scatterings, apart from 
the energy density e also the total particle number den¬ 
sity n is conserved, which are given by e = / p ^ P / P and 
n = f f p with the notation f p = f d d p/(2i:) d . The fact 
that they are conserved may be described by a continuity 
equation in momentum space, such as 


d 

ffj- ( W p/ P ) + V p ' j p = 0 


(Al) 


for energy conservation. Similarly, particle number con¬ 
servation is described by formally replacing w p —>■ 1 in 
the above equation and a corresponding substitution of 
the flux density. For the isotropic situation we can con¬ 
sider the energy flux A(k) through a momentum sphere 
of radius k. Then only the radial component of the flux 
density j p is nonvanishing and 


The first term in brackets comes from the scaling of 
the integral measure, the second from energy-momentum 
conservation for two-to-two scattering and the third from 
the three factors of the distribution function appearing 
in ((42]). 

Similar to the discussion in section |IVB| to be more 
general we may write for the scaling behavior of a generic 
collision term for (-vertex scattering processes 


C'(p) = |pf‘ C(l) 


(A7) 


in terms of the scaling exponent /q. One obtains along 
these lines 


jxi = (l — 2 )d — (Z + 1) — (i — 1 )k ■ (A 8 ) 

For the scaling properties of the energy flux we can then 
write 

AW °- 2WW2 + 1) /^IPl" S '".C(l).(A9) 

If the exponent in the integrand is nonvanishing, the in¬ 
tegral gives 

h.d+1+p.i 

m ~ (A10) 

d +1 + m 

Thus, scale invariance may be obtained for 

d+l + iH = 0. (All) 


/ Vp-jp = [ jp-dAp = (27 j) d A(k). 
J p J dk 


(A2) 


Since in this approximation uj p is constant in time, we 
can write with the help of (Al) and the kinetic equation: 


m = 


i 


2 d Tr d / 2 T(d/2 + 1 ) 


dp|p| 


>C 2 ~ 2 (p). 

(A3) 


This gives the scaling exponent for the perturbative 

l 


relativistic energy cascade: k = d — 


l- 1 


(A12) 


One observes that stationary turbulence requires in this 
case the existence of the limit 


lim 


C( 1) 


d+l+/ii —>-0 d -f- 1 jli 


= const 7 ^ 0 , (A13) 
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such that the collision integral must have a corresponding 
zero of first degree. Similarly, starting from the continu¬ 
ity equation for particle number one can study stationary 
turbulence associated to particle number conservation. 
This leads to the perturbative 


relativistic particle cascade: n 


= d — 


l + 1 

r ■ 


(A14) 


Accordingly, for quartic self-interactions we get k = 
d— 4/3 for the energy cascade and k = d— 5/3 for the par¬ 
ticle cascade. In the presence of a coherent field, when 
a 3-vertex can become relevant, the associated scaling 
exponents are k = d — 3/2 for the energy cascade and 
k = d — 2 for the particle cascade. For d = 3 this re¬ 
sult from stationary turbulence compares reasonably well 
with the approximate power law observed in section pll B| 
for the case of a self-similar evolution. For the corre¬ 
sponding perturbative nonrelativistic results we refer to 
the literature mmm- 


2. Nonperturbative regime: Strong turbulence 


The above perturbative description is expected to be¬ 
come invalid at sufficiently low momenta. In particular, 
the occupation numbers / p ~ jp| _K for k > 0 would grow 
nonperturbatively large in the infrared such that the ap¬ 
proximation (40) becomes questionable. To understand 


where the picture of weak wave turbulence breaks down 
and to compute the properties of the infrared regime, we 
have to consider nonperturbative approximations. For 
this purpose, we employ the vertex-resummed theory as 
in section [V] based on the expansion of the 2PI effective 
action in the number of field components to NLO. 

For a relativistic theory with dispersion ui p = 


i/p 2 + to 2 , the scaling assumption (A4) should be valid 


for sufficiently high momenta |p| m such that the dis¬ 
persion is approximately linear with w p ~ |p|. However, 
if a mass gap exists the infrared modes behave effectively 
nonrelativistic as explained in section |TV C[ In the follow¬ 
ing, we will analyze the two cases of a relativistic theory 
without mass gap and a nonrelativistic theory separately 
for comparison. 

We first look for relativistic scaling solutions in the 
absence of a mass gap. To be able to cope with occu¬ 
pancies of order ~ 1/A, we replace the perturbative colli¬ 
sion term by the vertex-resummed expression according 
to (87) with 


dp 0 dl°dq°dr 0 


I dn^°(p,i, q ,r) ^ ± (2ir)4 _ (J+1) 

x / S(p + l-q-r)p p pip q p r Xl s (p,l,q 1 r),(Al5) 

J lor 


where Ag ff is defined in ( 88 ). We emphasize that this still 
involves integration over frequencies as well as spatial 
momenta since no free-held form for the spectral function 


p p = p(p°,p) is used so far, and f p = f(p°, p). A cru¬ 
cial difference to the perturbative kinetic equation is the 
appearance of the vertex function u e ff(p°,p), which en¬ 
codes the emergence of a momentum-dependent effective 
coupling from the NLO corrections of the 1 /N expansion. 

We characterize the scaling properties of the spectral 
function p(p°, p) as in (67), which takes into account a 
possible anomalous dimension 77 . Accordingly, the scaling 
behavior of the statistical correlation function 


F(p°, p) = s 2+Kb F{s z p°,sp) 


(A16) 


will be described using a scaling exponent k s . This trans¬ 
lates with the definition for f p 1 into 

f(p°,p) = s^f(s z p°,sp). (A17) 


Using these definitions, one can infer the scaling behavior 
of v e g(p°, p), which is given in terms of t he “on e-loop” re¬ 
tarded self-energy n#(p°, p) as in section VB from which 
follows 


n«(p°,p) = s A U R (s z p°,sp) (A18) 


with 


A = 4— d — z + k s — 77 . (A19) 

If A > 0 one finds the infrared scaling behavior 

«eff(P°,P) = s~ 2A v e ff(s z p°, sp). (A20) 


(For A < 0 the effective coupling would become trivial 
with Deff ~ 1, on which we comment below.) Employing 
these scaling properties, (A15) gives 


J dn NLO (p,l,q,r) = s - 2 /C-*- 2 r, 
x J dfl NLC> (s z p°, s z l°, s z q°, s z r°; sp, si, sq, sr).(A 21 ) 


Following the procedure of section[Aj for any conserved 
quantity we can compute the flux through a momentum 
sphere k. Stationary turbulence solutions then require 
tha t th e respective integral does not depend on k. Similar 
to (A3), the flux for this effective particle number reads 


Mk) = ~ 2 -W(<i /2 + l) / ‘ , »l p l“' 1 C ' NLO tJ>>- 

(A22) 

The momentum integral can be evaluated along similar 
lines as before using the above scaling properties such 
that 


h,d— K s +z— 77 

A(k) ~ -A- C ,NLO (l) . (A23) 

w d- K s + z-r) v ’ y ' 

Therefore, scale invariance may be obtained for the Hug 

particle cascade: k s = d + z — 77 (A24) 



















19 


in the nonperturbative low-momentum regime. Similarly, 
for the scaling solution associated to energy conservation 
one finds, taking into account the additional power of p° 
in the respective integrand, the exponent for the mm 


energy cascade: k s = d + 2z — rj. 


(A25) 


With these solutions, we can now reconsider the above 


assumption that A > 0 by plugging (A24) or (A25) into 


(A19). Indeed, it is fulfilled under the sufficient condition 
that r) < 2. Taking into account that the anomalous 
dimension for scalar field theory is expected to be small 
for not too low dimension, we find k s ~ d + 1 for the 
particle cascade and k s ~ d + 2 for the energy cascade 
solution. 

The scaling of the effective occupation number distri¬ 
bution /(p), which depends only on spatial momentum, 
can finally be obtained from 


f (p) + l=J OO ^^ 0 p(P°,P) 


f(P°, P) + l 


^ 2W,P), 


(A26) 


using that f 0 dp 0 / ( 2 tt)p° p(p°, p) = 1 / 2 , which corre¬ 
sponds to the field commutation r elatio n in Fourier space. 
Of course, the scaling behavior (A16) for F(p°, p) may 


only be observed in a momentum regime with sufficiently 
high occupancies /(p) as implied in the above derivation. 
Then (A26) yields: 


/(P) = 


_ ki s +2— 2z 


f(sp) 


I — (d+2—z—rj) 


rel. particle cascade 
| — (d+2—97) re ^ ener gy cascade 


(A27) 


These estimates show that vertex corrections can lead to 
a strongly modified infrared scaling behavior as compared 
to the perturbative treatment of section EU 

We now turn to the nonrelativistic limit, which is out¬ 


lined in section IV and we present here the relevant 
changes as compared to the relativistic case. To this 
end, we consider a nonrelativistic iV-component complex 
scalar field theory and perform again the 1 /N expansion 
to NLO. Since already the relativistic scaling exponents 
(A27) indicate no explicit dependence on N at NLO - it 


can only enter indirectly via rj and 2 - this seems to be 
a very good starting point to understand also the single 
complex held case of Gross-Pitaevskii. 

Therefore, proceeding in the same way as for the rela¬ 
tivistic theory, with the corresponding scaling ansatz for 
the spectral and statistical two-point functions, leads to 


the very same solutions (A241 and (A25). A crucial dif¬ 


ference arises when the occupation number distribution 
is determined. Using the nonrelativistic definition for the 
distribution function, we have with the notation (701 in 
the absence of a condensate: 

1 Z 100 dn° 

/(p) + 2 =J -^F aa {p°,p). (A28) 


Comparison with the relativistic case shows that a differ¬ 
ence in the scaling behavior is cau sed b y the additional 
factor of ~ p° in the integrand of (A26 1 . Therefore, we 
find that 


/(P) = s K ‘ +2 ~ z f( S p) 

) |p| — ( d + 2 —^ 7 ) nonrel. particle cascade 
1 |p|V d+2+z-? d nonrel. energy cascade 


(A29) 


scales with one “z” difference than the relativistic solu¬ 


tion (A27l. This can have important consequences, such 


as the fact that the scaling exponent for the nonrela¬ 
tivistic particle cascade is now independent of the dy¬ 
namic exponent z describing the dispersion w p ~ |p|^. 
For d = 3 and 77 = 0 one has the scaling ~ l/|p| 5 . This 
result from stationary turbulence is not far from the ap¬ 
proximate power law ~ l/|p| K> described in section IB 
for the case of a self-similar evolution. 


Appendix B: Transport equations from quantum 
field theory 

In this section, we discuss the most relevant approx¬ 
imations involved in the derivation of the kinetic equa¬ 
tions used in the main text. We refer to Refs. mmm 
[SU i, ‘5Tj for further details. 


1. Nonrelativistic transport equations 

We consider the time evolution of the statistical and 
spectral propagators F and p = ip defined in 0 and 
(69), which obey 


f x ° 

[i°lc d x° - M ac (x)\F cb (x, y)= dz T, p ac (x, z)F cb {z, y) 

J to 


Jto 


dz^(x,z)p cb (z,y), 
(Bl) 

r° 

[ial c d x o - M ac (x)\p cb {x,y) = / dz T, p c (x, z)p cb (z, y), 

JyO 

(B2) 


where 

= &ab 


V 2 

2m 


+ | (f cc (x,x) + ip c {x)ip*(x)^ 

+ 9 (F a b(x, x) + ipa{x)tpt {x)j (B3) 

and ip a {x) = (ip a {x)) is the macroscopic field, computed 
from complex nonrelativistic quantum fields i/ji = ip and 

4>2=fo- 

Without approximating the self-energies, the above 
equations would be exact for Gaussian initial conditions. 
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The approximate expressions for the statistical and spec¬ 
tral parts of the self-energy X p and T, F used in this paper 
are given in appendix |B 3| An equation for the macro¬ 
scopic field completes the set of equations of motion in 
the case of a nonzero expectation value, which will not be 
used in the following. It can be found in reference m- 
The derivation of the transport equation for F (for p) 
follows standard procedures. One starts by switching the 
variables x and y in (Bl) (in (B2|) and subtracting the 


new equation from the original one. We change to Wigner 
coordinates X p = (x p + y p )/ 2 and s p = x p — y p and 
Fourier transform with respect to the relative coordinate 
according to 


g(X,p) = f ds p e^ g(x+ S -,X- |) 


(B4) 


Then one performs a gradient expansion, i.e. an expan¬ 
sion in orders of dx* and d p , keeping only the first order 
in this expansion. In particular, this implies that 

I d d+1 s xy e ip, * s *y I d d+l z f (x, z) g(z, y) 

* f(X,p)g(X,p) (B5) 

for any two functions /(x, y) and g(x , y). Writing A' 0 = t 
and dropping the X-dependence due to spatial homo¬ 
geneity, one arrives at the transport equations to leading 
order gradient expansion 

dF ab 
dt 

d Pab 


(t,p) = - i°lc [Fcd{t,p)^ p db {t,p) - p cd (t,p)i:% b {t,p)] 

(B6) 


dt 


(t,p)= 0. (B7) 


2. Relativistic transport equations 

To derive the transport equations for a relativistic 
scalar n-component theory, we start with the 2PI evo¬ 
lution equations for the statistical and spectral propaga¬ 
tors 

r x ° 

[S ac O x + Ml c (x)]F cb {x, y) = - dz E£ c (x, z)F cb (z, y) 

d to 
rv° 

+ / dzY, F c (x,z)p cb (z,y), 
J tQ 

(B8) 

r° 

\6 ac O x + M% c (x)]p cb (x, y)=- dz X p c (x, z)p cb {z, y), 

JyO 

(B9) 

where <p a (x) = (<p a (x)) and we defined the mass function 
A 


M ab(x) =S ab 


m 2 + — (f cc (x, x) + </> c (x)(/> c (x)) 


+ ^ (Fa b (x, x) + <j) a (x)<l> b (x)') . (BIO) 


The relativistic statistical and spectral functions F and p 
are defined in analogy to equations © and ( [69] ) by sub¬ 
stituting the operators t /> Q by scalar field operators (p a . 
Performing a gradient expansion to lowest order as ex¬ 
plained in the previous subsection, one obtains the trans¬ 
port equations 


= ~i [Z p c (t,p)F cb (t,p) - T, F c (t,p)p cb (t,p)] 
2p°^(t,p) = 0. (B12) 


(BH) 


We note that the spectral function p is time independent 
for both relativistic and nonrelativistic theories at this 
order of the gradient expansion. 


3. 2PI 1/A expansion to NLO 

We consider here the 2PI 1/iV expansion to NLO, 
where we consider the case of a vanishing macroscopic 
field, i.e. (f> a = 0 [35]. The corresponding equations for a 
nonzero field expectation value can be found in (37l . See 
reference [45j for an introductory presentation. 

At NLO the self-energies are given by 

^ab( x ) = B (l F {x,y)F ab {x,y) - ^ I p (x,y)p ab {x,y)\ , 


(B13) 


Kb( x ) = B ( IP ( X ’ y)Fab{x , y) + I F {x , y)pab(x, y) I, 


(B14) 

where B nr = —2g and B re \ = —A/3 N are the nonrela¬ 
tivistic and relativistic values for the prefactor B , 

r x ° 

! F (x, y) = n F (x, y)- dz I p (x, z)Ii F {z, y) 

J £q 
ry° 

+ / dzI F (x,z)H p (z,y), (B15) 

Jto 

r° 

I p {x,y)=Il p (x,y)- dzI p (x,z)U p (z,y), (B16) 
J v ° 


and 


B, 


n (x,y) = - — Tr 


F{x, y)F(y, x) + -p(x, y)p{y, x) 


n p (x, y)= - B Tr \p{x, y)F{y , x)]. 


(B17) 

(B18) 


Here F and p are to be understood as matrices. Next, the 
initial time is sent to the remote past to —> — oo, and we 
introduce retarded and advanced quantities as I R (x, y) = 
0(x° — y°)I p (x,y), I A = —6 (y° — x°)I p (x,y) and similar 
for and n l . Then (B16) is Fourier transformed with 
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respect to the relative Wigner coordinate to obtain at 
lowest order in the gradient expansion 


I R (X,p) 


n R (x, P ) 
i + n fl (x,p)' 


(B19) 


Using I p = I R — I A and proceeding in the same way with 
(B15), one obtains 


with 


I F (X,p)=v eS (X,p)n F (X,p) (B20) 

I p (X,p) =v eS (X,p)n p (X,p) (B21) 


V eS (X >P )- | 1 + n fl (X)P )|2 ( B22 ) 


and 


n R (X,p)= - B [ F ab (X : q-p)G R (q). 

J a 


(B23) 


Inserting these expressions into the transport equation 
(B 61 and (Bll), one obtains the evolution equations (72) 


and (871, respectively. 


The effective coupling can be further simplified in the 
approximation with (84), the identity 


C£M = lim j F 


dio -ip ba (u, p) 


oj — — ze 


(B24) 


and the on-shell spectral functions (83) and (66) for the 


nonrelativistic and relativistic theories, respectively. In 
this way, one arrives at the quasi-particle expressions (58) 
and ([55|. 


Appendix C: Extracting scaling exponents for a 
self-similar evolution 


In sections El and nn we observe that the nonrelativis¬ 
tic and relativistic N = 2 and N = 4 theories evolve in 
a self-similar way, characterized by ([l]). We extract the 
pair of scaling exponents (a, /3) and estimate their uncer¬ 
tainties from this self-similar evolution. For our analysis, 
we use a similar least-squares method as introduced in 
reference [T3] . In the following, we describe the method 
used in the present paper. 

At first, we quantify the deviation from the self-similar 
evolution of the distribution function. Using the notation 


p = |p|, for each pair (a,f3) the spectra at several times 
are rescaled according to 

fresc{t,p) = (t/t Iei )~ a f(t, {t/treiypp) . (Cl) 

We use the distribution at the earliest time t re f as a 
reference and the distributions at the N com later times for 
comparison. A perfectly self-similar evolution ([I]) implies 
Xf(t,p) = f r esc(t,P)—f(tref,P) = 0 such that the rescaled 
distribution function becomes time independent. 

However, this equality is in general only true for the 
correct set of scaling exponents and is in practice violated 
due to statistical uncertainties of the data and systematic 
deviations from the perfect scaling behavior. Minimizing 
the deviations yields the best fit for the scaling expo¬ 
nents, and the study of the distribution of the deviations 
gives an estimate for the uncertainties of the exponents. 
We quantify the deviations by 

2 , gs 1 f d(log(p)) (Af(t k ,p)/f(t re{ ,p)) 2 
X [a,P ’ N com 2^ f d(log(p)) 

(C2) 

where we sum over all relative deviations for each com¬ 
parison spectrum. We use integration over d(log(p)) to 
increase the sensitivity at low momenta. The upper limit 
of integration is given by the highest momentum included 
in the inverse particle cascade, i.e. the last point of the 
approximate 1 /p K> power law of the reference distribu¬ 
tion function. Since the spectra are binned in momen¬ 
tum space, the integrals translate to sums over momenta 
f d(log(p)) —¥ Ty^logipi+i/Pi), where p i+ 1 > p t are 
the discrete momenta of each testing spectrum k with 
possibly a different number of bins n*,. To compute the 
difference A/(t/-,pi), we linearly interpolate momenta of 
the reference spectrum to coincide with the discrete mo¬ 
menta of the rescaled spectra. 

The deviation x' 2 ( a - 0) is minimal for the best fit of the 
scaling exponents a and /3. In analogy to reference 
we define the likelihood function of a given set of expo¬ 
nents (a, (3) as 


W(a,0) 


AT 6Xp 


2 x 2 (d,/3)J ’ 


(C3) 


where M is a normalization constant such that the inte¬ 
gral f da d/3W = 1. Integrating W{a,j3) over one of the 
exponents provides a marginal likelihood function for the 
other one, e.g. W(a ) = / d/3 W(a, (3). Approximating 
the marginal likelihood functions with Gaussian distribu¬ 
tions, we obtain an estimate for the standard deviations 
cr a and ap, while the means are still given by a and (3. 
Our thus measured scaling exponents are finally written 
in the form 


a = a±crQ,, f3 = [3 ± ap . (C4) 
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